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Abstract. 

We present results from our ongoing follow-up observations of double white dwarf binaries detected in the ESO SN la 
Progenitor SurveY (SPY). We discuss our observing strategy and data analysis and present the orbital solutions of five close 
double white dwarf binaries: HE0320-1917, HE1511-0448, WD0326-273, WD1013-010 and WD1210+140. Their peri- 
ods range from 0.44 to 3.22 days. In none of these systems we find any spectral lines originating from the companion. This 
rules out main sequence companions and indicates that the companion white dwarfs are significantly older and cooler than 
the bright component. Infrared photometry suggests the presence of a cool, helium-rich white dwarf companion in the binary 
WD 0326-273. We briefly discuss the consequences of our findings for our understanding of the formation and evolution of 
double white dwarfs. 
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1. Introduction 

Most binaries in which the two components are close enough 
together to interact at some stage, will en d their life as close 
double white dwarfs (e.g. Webbink 1984). They will experi- 
ence at least two phases of mass transfer. In at least one of 
these phases, a so-called common envelope will be formed, in 
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which the first formed white dwarf spirals into the envelope of 
its companion (e.g. lNelemans et alJ200lh . 

Although this population of close double white dwarfs was 
expected on theoretical grounds (Web binkll98l. it was only i n 
1988 that the first system was discovered JSaffer et al1ll988t) . 
In t he 1990s the number of known s ystems increased to 18 
(see lMarsrJ2000llMaxted et all2000al) . Detailed knowledge of 
this population is important for a number of reasons. The ob- 
served double white dwarfs carry important information about 
the preceding binary evolution; in particular the mass ratio of 
the components reflects the change in orbital separation in the 
mass transfer phases and puts constraints on the physics and 
applic ability of the common envelope phase ( Nele mans et alJ 
2000). Furthermore double white dwarfs are so abundant that 
they form an important background for gravitational wave ra- 
diation experiments l ike LISA, as well as be ing guaranteed de- 
tectable sources (e.g. lNelemans et alJl2004 . Last but not least 
double white dwarfs play an important role in the still ongoing 



2 



G. Nelemans et al.: Binaries discovered by the SPY project 



Table 1. Overview of the follow-up observations, giving dates, telescopes, instruments and observers. For details of the set-ups 
of the different instruments, see text. 
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debate on, or rather quest for, the nature of the progenitors of 
type la supernovae. 

Type la supernovae (SNe la) are among the best stan- 
dard candles to d etermine the cosmological parameters (see 
iLeibundeutl 1200 ll) . Yet t he nature of the SNe la progeni- 
tors remains unclear ("see iLiviol 1200 ll for a review). One of 
the p rogenito r models is the merger of two white dwarfs 
Jlben & Tutukovll984 . Some explosion calculations show that 
merging doub le white dwar fs do not lead to an explosive event 
(e.g. ISaio & Nomotolll998h even when rotation is taken into 
account ( Sai o & Nomotol 120041) . w hile other calculation s in- 
cluding rotation do get explosions JPiersanti et allEooil) . The 
current explosion models are still too simplistic to rule out 
double white dwarfs as progenitors. For the use of SNe la as 
standard c andles the nature of the progenitors is crucial (e.g. 
iHoeflich et al.ll 19981) . One basic question that needs to be an- 
swered for any possible progenitor scenario is whether there 
are enough of them. 

To give a definite answer to the question whether there 
exist enough merging double white dwarf with a mass above 
the Chandrasekhar limit we embarked on a spectroscopic sur- 
vey of white dwarfs with the UVES spectrograph at the ESO- 
VLT (ESO SNIa Progenitor Survey, INapiwotzki etaDl2001ah 
and have observed more than 1000 white dwarfs (which corre- 
sponds to about 75 per cent of the known white dwarfs acces- 
sible by the VLT with magnitudes brighter than B = 16.5). For 
each object we took two spectra to check for binarity, in order 
to provide good statistics of the double white dwarf population, 
including SN la progenitor candidates. A total of 875 stars were 
checked for radial velocity variations. SPY detected « 100 new 



double white dwarfs (see INapiwotzki et al.1 2004b). There are 
16 double-lined systems and 18 non-DA white dwarf binaries. 

The sample furthermore turned out to be contaminated with 

j 1 

some sdB stars, 31 of which are binaries (Napiwotzki et al. 
2004a). The SPY project obtains 2 to 3 spectra per object, 
enough to determine its binarity and pick out the 1 in 10 or 
so white dwarfs that reside in a close binary. However, in or- 
der to obtain orbital solutions of all the binaries discovered in 
the SPY project, we started a systematic study of all new bina- 
ries, mai nly with smaller telescopes. Earli e r results were pub - 
lished in INapiwotzki et all feOOlbl l2002bl) : Ikarl et all J2003I) . 
An overview of the follow-up observations up to Dec 2004 is 
given in table [0 In this paper we report results for the first 
five single-lined double white dwarfs for which we found the 
orbital solutions. Results on other types of binaries (double- 
lined systems and sdB binaries) will be presented in separate 
papers. In Sect. [2] we discuss our observations, data reduction 
and analysis, in Sect. [3] we present the orbital solutions for the 
five systems and in Sect.[5]we briefly compare our findings with 
population synthesis calculations. 

2. Method 

2.1. Observations and data reduction 

Data were taken with a number of different tele- 
scope/instrument combinations (see Table 0. The details 
of the different setups are as follows. 

UT2/UVES: data taken with the UVES echelle spectrograph 
on UT2 of the 8.2m Very Large Telescope of the European 
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Southern Observatory at Paranal with the following setup: 
dichroic 1, central wavelengths 3900 and 5640 A giving an al- 
most full spectral coverage from 3200 to 6650A with two small 
(~ 80A) gaps around 4580 and 5640A. Slit widths varied be- 
tween 1" and 1.5" for the follow-up observations and 2.1" for 
the survey spect ra, yielding a spectral re solution of 0.36A at 
Ha or better (see Napi wotzki et all2001al) . 

3.5m/TWIN: data taken with the double arm TWIN spectro- 
graph on the 3.5m telescope at Calar Alto. Only data from the 
red arm are used in this paper. The setup consisted of the 230 
mm camera with the T06 grating (1200 grooves mirT 1 ) and a 
SITe-CCD giving a dispersion of 0.55 A/pix. 

INT/IDS: data taken with the Intermediate Dispersion 
Spectrograph in the 2.5m Isaac Newton Telescope on La 
Palma. The 500 mm camera with the 1200R grating centred on 
Ha was used, giving a dispersion of 0.39A/pix. 

WHT/ISIS: data taken with the double arm ISIS spectrograph 
on the 4.2m William Hershel Telescope on La Palma. Only data 
from the red arm are used in this paper. The setup was the 
1200R grating centred on Ha, together with the MARCONI 
CCD giving a dispersion of 0.23 A/pix. 

NTT/EMMI: data taken with the EMMI spectrograph on the 
New Technology Telescope at La Silla in the red medium dis- 
persion mode (REMD). We used grating #6, giving a dispersion 
of 0.2 A/pix on the MLT/LL CCD. 

Data reduction was done using standard IRAF or MIDAS 
tasks. Bias subtraction was done using the overscan region, 
while flatfielding was done using daytime dome flats. Spectra 
were extracted using optimal extraction, except for the survey 
spectra. Wavelength calibration is based on CuArNe arcs taken 
at the same position as the object before or after each exposure 
(for the IDS/ISIS data) and ArTh for the TWIN data. Details 
of the UVES data reduction will be described in a forthcoming 
paper on the SPY radial velocity catalogue (Napiwotzki et al, 
in prep). 

2.2. Radial velocity measurement and period 
determination 

Based on previous experience we have used two different tech- 
niques of measuring the ra dial velocities of the white dwarfs. 

The first is described in Napiwotzki et al . (2002b) and uses 
functions within the MIDAS context to fit a combination of a 
Gaussian plus a Lorentzian to the spectra. The radial veloc- 
ities are calculated from the central wavelength of the fitted 
Gaussians. The velocities and observation times were later cor- 
rected to heliocentric values. In order to determine the orbital 
parameters radial velocity curves for a range of periods are fit- 
ted to produce a "power spectrum" indicating the fit quality as 
function of period. 

We also used the technique developed by Marsh and collab- 
orators, which utilizes a multi-Ga ussian fitting technique to fit 
the shape of the spectral lines (see M arsh et alJll995h . A fixed 
template consisting of a number of different Gaussians (typi- 
cally two or three), is fitted to the individual spectra in order 



to measure the radial velocity. A floating mean period search is 
used to find the orbital solution by minimizing the x 2 of the fit 
of a sinusoida l curve to the radial velocity measurements (see 
iMorales-Rueda et"ai1l20 05. for a discussion of this technique). 
This yields the best fit period, radial velocity amplitude, veloc- 
ity offset (in principle the system velocity, or y-velocity, but not 
corrected for gravitational redshift), the zero-point (where the 
star moves from the blue to the red) and the mass function. 

In order to cross check the two analysis methods we com- 
pared the resulting radial velocity measurements in detail for 
two binaries (HS1 102+0934 and HS1 102+0032, details of 
which will be published later) and found good agreement. Most 
values agree within 1 sigma, the maximum deviation was 1 .44 
sigma. All radial velocities presented in this paper are obtained 
with the second method. 

Before discussing the individual objects we spend a few 
words on the error analysis. Following earlier practice with ra- 
dial velocity data of double white dwarfs, we assume there are 
systematic errors of a few km s _1 that are due to imperfections 
in the setup. We therefore added 2 km s _1 quadratically to the 
errors that we obtain by propagating the errors on the spectra. 
We found that there were some additional errors that arise from 
combining data from different telescopes/instruments. We did 
not include any additional error term to account for this, but 
we chose to keep the errors as they are, which in some cases 
results in rather large values of^f 2 . The errors on the derived 
quantities might thus be a bit optimistic. However, none of the 
derived periods is uncertain due to this effect. 

3. Results 

Table|5]summarizes the results for this first set of double white 
dwarfs. Listed are the orbital period (P rt>X the radial veloc- 
ity semi-amplitude (K), the system-velocity (y), the zero-point 
(To, where the star moves from the blue to the red), the x 2 of 
the fit, the number of data points (N), the difference of the x 2 
with the^- 2 of the second best fit period (Ay 2 ), the mass of the 
visible white dwarf (Mi) and the minimum mass of the com- 
panion white dwarf, from the mass function, M\ and assuming 
an inclination of 90 degrees (M2 y m m )- 

The masses M\ of the visible component were determined 
from a model atmosphere analysis of th UVES survey spec- 
tra (Voss et al. 2005, in prep.). In short, we determined tem- 
perature and surface gravity by fitting model spectra to the 
observed hydrogen Balmer lines of the DA white dwarfs. 
Depending on the temperature of the star we applied LTE or 
NLTE model atmospheres (cf. iKoester et alJl200ll for a de- 
scription ). Finally we interpola t ed in the white dwarf cooling 
tracks of iBenvenuto & Althausl i 19991) to determine the mass. 
The results for the systems are discussed below individually. 

3.1. HE0320-1917 

HE0320-1917 is relatively cool white dwarf of about 12,000 
K, with a mass of 0.29 M . It was selected from the 
Hamburg/ESO survey (see ichristlie b et alJl200ll) . The peri- 
odogram of HE0320-1917 is shown in Fig.^ The best fit pe- 
riod is about 0.865 d (see Table 0. The radial velocity mea- 
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Table 2. Overview of results of the follow-up observations. The columns give orbital period (P 0! b), the radial velocity semi- 
amplitude (K), the gamma-velocity (y), the zero-point {To, where the star moves from the blue to the red), the^ 2 of the fit, the 
number of data points (N), the difference of the;^ 2 with the^- 2 of the second best fit period (Ay 2 ), the mass of the visible white 
dwarf (M\ ; from Voss et al. in prep.) and the minimum mass of the companion white dwarf, from the mass function, M\ and 
assuming an inclination of 90 degrees (M2 itn j n ). 
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Fig. 1. Top: periodogram for HE0320-1917 showing the^ 2 as 
a function of the orbital period in cycles per day. Due to the 
availability of measurements over a long time span the x 2 min- 
ima are so narrow that they appear as vertical lines. Bottom: the 
Radial velocity data folded on the best fit period. 



surements folded on the best fit period are shown in Fig.^ The 
minimum mass of the unseen companion is 0.35 M . 



Fig. 2. Periodogram for HE15 11-0448 in cycles per day (top) 
and data folded on the best fit period. 



3.2. HE1 511 -0448 

HE15 1 1 -0448 was again selected from the Hamburg/ESO sur- 
vey and is a very hot white dwarfs (50,000 K), with a mass of 
0.48 M . The line profile shows an emission core, so this ob- 
ject was fitted with three Gaussians, one for the emission core. 
The periodogram of HE151 1-0448 is shown in Fig. 13 The or- 
bital period is 3.22 days. There are two other solutions for the 
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Fig. 3. Periodog 
and data folded 



ram for WD0326-273 in cycles per day (top) 
on the best fit period. 



-25 - 



-25.5 



-26 



v 



-25 



-26 



single DA 



i i i i I i 
composite DA+DC 



DA 



25.5 - DC 



5000 



10" 



1.5xl0 4 
A/A 



2x10* 



Fig. 4. Top: Model flux of a 9300 K DA white dwarf compared 
to the observed fluxes (converted from magnitu des). For clarity 
we show only the JHK measurements of Be rgeron et alJj200ll) . 
Bottom: Comparison of the best fit composite DA+DC spec- 
trum to the observed fluxes. The individual contributions of the 
DA and DC components are shown as dashed lines. 



period (0.8 and 1.3 d), but they are very unlikely. The folded 
radial velocity curve is shown in Fig. [2] 

3.3. WD0326-273 

WD0326-273 (GJ 1060A, L587-77A) was discovered by 
Luyten as a high proper motion white dwarf with an M(3.5) 
companion at ~7 arcsec (e.g. iLu vtenll 19491) . The velocity dif- 
ference of the two stars has been used to determine the grav- 
itational redshift of the white dw arf (Koester 1987), yielding 
a mass of 0.65 M^. IWegner & Reidl dl99ll) on the other hand 
measured a blue shift and discarded that system as not phys- 
ically associated. However, both these results are, of course, 
influenced by the orbital velocity. 

Spectral fits to the SPY UVES spec tra yield a ma ss of 
0.46M o (Voss et al. 2005, in prep.). lBergeron et alJ J200ll) find 
a mass of 0.35 M based on photometric and parallax data, 
but a higher mass using spectroscopic mass determination. The 
analysis of the photometric data was flawed, because Bergeron 
et al. adopted V = 14.0, a value which can be traced back to a 
photographic measurement bv lLuvtenl ill 9571) . A more accurate 
value (y/V=13.56) is available from photoelectric observations 



in the Stromgren system bv lBessell & Wickr amasinghe] lll978l) . 
We will discuss evidence for a remaining infrared excess below. 

Gravitational redshift: With our orbital solution we are now 
in a position to determine the real gravitational redshift of 
the visible DA white dwarf. Due to a (fortunate) target con- 
fusion during the SPY survey observations one UVES spec- 
trum of the M star companion L587-77B was taken. The 
Balmer lines are in emission and we measured the radial ve- 
locity from a simple Gaussian fits of the Ha and HB emis- 
sions. The measured value VMstar = 47 + lkms~' translates 
into a redshif t of 23kms~ 1 . Apply i ng the white dwarf cool- 
ing tracks of iBenvenuto & Alfhausl ([T999) we derive a white 
mass of 0.5 1M , slightly higher than the spectroscopic result. 
Let us note that taken at face value our radia l velocity deter- 
minati on of the M dwarf is discrepant to the IWegner &~ Reid 
d!99lh measurement of VMstar = 17 ± lOkms However, the 
authors caution that the errors of their absolute radial velocity 
measurements are probably larger than the purely instrumental 
error. 

The cool white dwarf companion: The quality of the model 
atmosphere fit of WD 0326-273 is inferior to the quality typ- 
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Fig. 5. Periodogram for WD1013-010 in cycles per day and Fig. 6. Periodogram for WD1210+140 in cycles per day and 



data folded on the best fit period. 



data folded on the best fit period. 



ically achievable for white dwarfs in this parameter range 
(Teff = 9300K, Voss et al. 2005, in prep.). Closer inspec- 
tion of the optical and IR photometry revealed that a signif- 
icant red/infrared excess is present (Fig. @] top). A plausible 
explanation is a flux contribution from the apparently invisi- 
ble cool white dwarf companion of the brighter DA. We col- 
lected the availab le photometric measurements (Stromgren b 
and y from iBessell & W ickramasinghel dl978l) : J,H,K from 
iBer geron et alJ feOOll) : / and ./ from DENIS: J,H,K from 
2MASS) and tried to fit these with a composite spectrum. 
We fixed the temperature of the bright DA component at the 
spectroscopic value of 9,300 K and adopted equal masses for 
both white dwarfs. Synthetic absolute magnitudes and colours 
for pure hydrogen a nd pur e helium atmospheres were taken 
from Bergero n et alJ ([1995). Fluxes were co-added and nor- 
malised to the measured V flux. A x 1 analysis yielded a best 
fit for a helium-composed companion with a temperature T e ff = 
6 150±650K (displayed in the bottom panel of Fig. |4}. This cor- 
responds to a featureless DC white dwarf. Fitting the colours 
for the hydrogen-rich models yielded a slightly worse fit at 
T e « = 7300 + 350 K. However, this solution can be ruled out, 



because it predicts a relatively strong Ha line for the secondary, 
which is clearly not present in the observed spectra (see Fig.0. 

The inner binary (assuming the M3.5 star is physically as- 
sociated in a very long period outer binary) has an orbital pe- 
riod of 1.88 d (Fig. There are no other periods possible, 
given the relatively large number of spectra we have of this 
source. The large amplitude of the radial velocity curve shows 
that the companion star also has to be rather massive (more 
massive than 0.59 M , see Table|2}. 

3.4. WD1013-010 

WD1013-010 is a cool white dwarf (8,000 K) with a mass 
of 0.44 M . The periodogram for WD 101 3-010 is plotted in 
Fig. |3 The best fit period is 0.437 d. The folded radial velocity 
curve is shown in Fig. [5] The next best period, with^- 2 =51 
has virtually the same period and thus overlaps with the best fit 
period in Fig.|5] 
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Table 3. Properties of the observed double white dwarfs 




Velocity (km/s) 

Fig. 7. Quadrature spectra of WD0326-273, after shifting out 
the best fit sinusoidal radial velocity curve and combining the 
spectra in the phase intervals 0.15-0.35 (top) and 0.65-0.85 
(bottom). Any signs of the companion would show up as an 
asymmetry in these two spectra. 



3.5. WD1210+140 

WD1210+140 is another hot white dwarf (32,000 K) of very 
low mass: 0.23 M Q . It is the lowest mass white dwarf known in 
a double white dwarf binary, and one of the lowest m ass white 
dwarfs known. Interestingly, Max ted & Marshl (ll999) observed 
WD1210+140 and, despite the large radial velocity amplitude 
of the system, concluded is is not a binary, but only based 
on two (consecutive) spectra. We included their two measure- 
ments (taken in 1998) in our fit. This results in a unique solution 
for the system, with orbital period of 0.64 d. (Tabled- The ra- 
dial velocity curve, fo lded on the best period is shown in Fig. [6] 
The measurements of iMaxted & Marshl i 199 9) are included as 
open triangles. The mass of WD1210+140 is virtually identi- 
cal to the mass of the sdB star in the 0.6 d binary HD 1881 12, 
in whi ch the c ompanion is a massive (M > O.73M ) compact 
object jHeber et alJ2003h . 



4. Is there any sign of the companions? 

In order to search for signs of weak absortion lines of the com- 
panions, that potentially could be used to determine the mass 
ratio, we have to inspect the spectra around quadrature (phases 
0.25 and 0.75) when the two spectra are most widely separated. 
In order to obtain high quality spectra to look for weak spectral 
features one generally combines spectra with similar phases. 
However, in general our single high singnal-to-noise ratio VLT 
survey spectra, that were used to measure the initial radial ve- 
locity offset are more suitable for this. In the survey spectra 
none of the systems discussed here showed any sign of absorp- 
tion lines of the companion. Only for WD0326-273 we can 
improve on these results as many of the follow-up spectra were 
taken with the VLT as well. We therefore combined the VLT 
spectra in the phase range 0.15-0.35 (4 spectra) and 0.65-0.85 
(9 spectra) after shifting out the radial velocities according to 
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0.35 


>0.23 
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0.48 


>0.46 


18 


WD1704+481 


0.14 
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>0.18 
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WD 1824+040 
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0.43 


>0.52 


8,17 


WD2032+188 


5.08 


0.41 


>0.47 


8,17 


HE2209-1444 


0.28 


0.58 


0.58 


16 


WD233 1+290 


0.17 


0.39 


>0.32 
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our best fit orbital solution. The resulting spectra are shown in 
Fig-El no additional absorption lines are seen, but see Sect. 13. 31 
for our arguments that there must be significant flux from a 
featureless companion. 

5. Discussion 

In total there are now 24 double white dwarfs known. The 
parameters of all systems are shown in Table [3] including 
the new systems discussed in this paper. Of the new sys- 
tems WD0326-273 and WD1210+140 are particularly inter- 
esting, having rather high (0.5 1 M ) and very low (0.23 M ) 
masses respectively. The difference in the approach between 
the earlier work (targeting low-mass white dwarfs) and the 
SPY project is starting to become clear. Of the four clear C/O 
white dwarfs (taken as M > 0.5M G ) three were found by 
SPY. There are a few changes compared to earlier compila- 
tions like this (e.g. IMaxted & Marshlfl99^: iNelemans & Toutl 
l2005t) . in particular some of the masses of the white dwarfs 
have been redetermined. There are als o some small changes 
in the masses given in new papers IMaxted et alJ |2002b; 
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log P(hr) 



Fig. 8. Period - mass (Mi) distribution for double white dwarfs. 
The systems presented in this paper are shown with the open 
circles. 



iMorales-Rueda et"afll2005l) . In the table we also list the min- 
imum companion star mass for the single-lined spectroscopic 
binaries, based on the mass function and assuming an inclina- 
tion of 90 degrees. We can now use the above data in order 
to compare with theoretical models of the local population of 
double white dwarfs to see if the new systems change the con- 
clusions of earlier comparisons. We compare the data with a 
slightly updat ed version of th e population synthesi s calcula- 
tion o f INelemans et alJ ll200ll) . as detailed in Nele mans et alJ 
Earlier population synthesis calculation s for double 
white dwarfs are reported e.g. in iTben et.alHl997fc : lHar] fl998). 

In Fig.[S]we plot the distribution in the period - white dwarf 
mass plane, for the brighter white dwarf in each binary (Mi). 
All observed double white dwarfs are plotted on top of the 
theoretical distribution for a magnitude limited sample with 
(V < 20). A few things should be noted before conclusions 
based on this diagram can be drawn. The most important is 
that the observed systems are not (yet) an unbiased sample. 
The SPY project will provide a much more unbiased sample, 
but only when the majority of the discovered binaries have 
known periods. More specific, the observed systems are still 
dominated by the low-mass white dwarfs that were targeted by 
Marsh and collaborators. 

Taking the plot as it is, two things can be concluded. The 
first is that the new observations occupy similar regions in the 
diagram and thus conform to the general picture that was seen 
before and which can be explained adequately by the popula- 
tion synthesis calculati ons. It also confirms the earlier finding 
( INelemans et aklEoOlh that the lowest mass white dwarfs do 
not keep a thick hydrogen envelo pe in which resid ual burn- 
ing continues, as was suggested bv lDriebe e t al. ( 1998), which 
would predict many systems with very low masses (although 
WD1210+140 shows that some systems with M < O.25M 
do exist). This might be due to thermal in stabilities in which 
the hydrogen is lost (e.g. Drie be et al Jl999h . However, we cur- 
rently use a cooling prescription based on the models o fH ansen 
( 1999), as detailed in Nelem ans et alJ J20041) . which assume a 



log P(hr) 



Fig. 9. Period - mass ratio (M1/M2) distribution for double 
white dwarfs. The systems presented in this paper are again 
shown with the open circles. The systems for which only a 
lower limit to the mass of the companion is know are shown 
with the arrows. They are drawn from the maximum mass ratio 
(for M2 = M2 jm in), to the point where there is a 75 per cent 
probability that the mass ratio lies above this point. The large 
circle is plotted at the position of the median inclination (60 
degrees). Open circles again are the systems presented in this 
paper. In the theoretical distribution Mi is selected to be the 
mass of the brightest object of the two (based on age and cool- 
ing, generally the second formed white dwarf). 



constant hydrogen envelope mass as function of the mass of the 
white dwarf. This prescription enhances the presence of sys- 
tems with mas ses above 0.5 Mn i n Fig. [8] as compared to the 
same plot in Nele mans et alJ J200 lh . more consistent with the 
observations. The second conclusion is that the apparent gap 
in the period distributio n betw een 0.5 and 1 day, which was 
noted in lNelemans et alJ J200 ll) indeed was due to small num- 
ber statistics and is filling in. 

In Table[5]we also included the minimum mass of the com- 
panion for the single-lined systems. The mass ratios of dou- 
ble white dwarfs 1 have led to the suggestion that the standard 
description of the common-envelope phase in close binaries 
needs revision (INelemans et all200d). The population synthe- 
sis calculations of lNeleTnans*et alJj200ll) are based on this sug- 
gested revision. A recent reassessment of the need for this revi- 
sion, based on the newly discovered double -lined double whit e 
dwarfs confirms the earlier findings ( Nelemans & Tout l2005t) . 
In these analyses only double-lined systems were taken into 
account. However, in Fig. [9] we compare the theoretical distri- 
bution in the period - mass ratio (Mi IM2) plane with all obser- 
vations, where for the single-lined systems we show a plausible 
range of mass ratios, based on the mass limit on the companion 
star. For each system we plot an arrow that starts at the max- 
imum mass ratio (for the minimum companion star mass, i.e. 
for 90 degrees inclination). The arrow extends to the mass ratio 



1 Traditionally taken to be the mass of the brighter ov er the mass of 
the dim mer component (Mi /M2), e.g. llben et al T <1997l) : lMa"xtedetalJ 
1 2000b) 
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for an inclination of 41 degrees, which corresponds to a 75 per 
cent probability that the actual mass ratio lies along the arrow. 
The large symbol is plotted at the mass ratio for the median 
for random orientations, i.e. an inclination of 60 degrees. In the 
theoretical distribution for each system the brightest of the two 
white dwarfs (in the V band) is determined, based on the age 
and the (mass dependent) cooling curve and its mass is taken 
as M\ . This generally is the second formed white dwarf. 

Although the general agreement between the observations 
and the theory is good, it clearly is not as good for the single- 
lined systems as for the double-lined systems (shown as sym- 
bols with error bars), for which the theory was developed. In 
particular the mass ratios for the longer period systems seem to 
be lower than predicted. Detailed investigation of this discrep- 
ancy is beyond the scope of this paper. 
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Appendix A: Radial velocity measurements 



Table A.l. Radial velocity measurements of WD0326-273 
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